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ABSTRACT 

Water is one of the key molecules in the physical and chemical evolution of star- and planet- 
forming regions. We here report the first spatially resolved observation of thermal emission of 
(an isotopologue of) water with the Plateau de Bure Interferometer toward the deeply embed- 
ded Class protostar NGC 1333-IRAS4B. The observations of the H* 8 3 U - 2 2 , transition at 
203.4 GHz resolve the emission of water toward this source with an extent of about 0.2" corre- 
sponding to the inner 25 AU (radius). The H* 8 emission reveals a tentative velocity gradient 
perpendicular to the extent of the protostellar outflow/jet probed by observations of CO rota- 
tional transitions and water masers. The line is narrow « 1 km s -1 (FWHM), significantly less 
than what would be expected for emission from an infalling envelope or accretion shock, but 
consistent with emission from a disk seen at a low inclination angle. The water column density 
inferred from these data suggests that the water emitting gas is a thin warm layer containing 
about 25 MEarth of material, 0.03% of the total disk mass traced by continuum observations. 

Subject headings: astrochemistry — stars: formation — planetary systems: protoplanetary 
disks — ISM: abundances — ISM: individual (NGC 1333-IRAS4B) 



1 Based on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported by INSU/CNRS (France), 
MPG (Germany) and IGN (Spain). 

2 Current address: Centre for Star and Planet Formation, Natural History Museum of Denmark, University of Copenhagen, 
0ster Voldgade 5-7, DK-1350 Copenhagen K, Denmark 
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1. Introduction 



Water is one of the most important molecules in star-forming regions: it is a dominant form of oxygen, 
is important in the energy balance, and is ultimately associated with the formation of planets and emergence 
of life. Thus, following the water "trail" from collapsing clouds to protoplanetary disks is a fundamental 
problem in astronomy and a strochemistry. In the c old and quiescent regions the gaseous water abundance is 
low, only 10~ 9 - 10~ 8 (e.g. jBergin & Snellll2002r) . but in regions with intense heating or acti ve shocks, its 
abund ance can reach 10~ 4 with respect to H2 - comparable to or higher than that of CO (e.g. jHarwit et al. 
1998|). Which mechanism is most important for regulating the H2O abundance in low-mass protostars is 
still heavily debated, ho wever. Is it passive heating of the collaps ing envelope by the accretion luminosity 
from forming stars (e.g., Ceccarelli et al. 1998k LMaret et al. 2002 ). or shocks caused either b y protostellar 
outflo ws (e.g., Nisini et al. 19991) or related to ongoing accretion onto circumstellar disks ( Watson et al. 
20071) ? H2O is also a key molecule in the chemistry in regions of star formation: in large parts of the cold 
and dense e nvelopes around low-mass protostars, H? Q is the dominant constituent of the icy mantles of dust 
grains (e.g.. lWhittet et al.lll988l : iBoogert et al.ll2008r) . Its evaporation at temperatures higher than 90-100 K 
determines at what point water itself and any complex organic molecules, formed in these ice mantles, are 
injected into the gas-phase. 

This discussion has received new fuel with the detection of surprisingly strong highly excited H2O 
lines at mid-infrared wavelengths with the Spitzer Space Telesc ope's infrared spec trograph, IRS, toward 
one deeply embedded Class protostar, NGC 1333-IRAS4B, by Iwatson et all (120071) . Based on the high 



critical density of the observed lines and temperature (~ 170 K). [Watson et al.l argue that the water emission 
observed toward this source has its origin in an accretion shock in its circumstellar disk. Those data could 
not spatially or spectrally resolve the water emission, however. 

In this letter we present observations at 203 GHz of the Hj, 8 isotopologue at high-angular resolution 
(0.5") using the Institut de Radioastronomie Millimetrique (IRAM) Plateau de Bure Interferometer to de- 
termine the origin of water emission in low-mass protostars. This isotopic line is a useful tracer of H2O as it 
can be detected and imaged at high angu l ar resolution from the gr ound under good weather conditions (e.g., 
Jacq et al.ll 19881 : iGensheimer et al.ll 19961 : Ivan der Tak et al.ll2006f) . Its upper level energy of 203.6 K is well 
matched to the observed excitation temperature of water seen by Spitzer. 



2. Observations 



We observed NGC 1333-IRAS4B (IRAS4B in the following; a=03 h 29 m 12!00, 5=+31 o 13'08'/l [J2000] 
IjQrgensen et al. 2007a ] 1 using the 6 element Institut de Radioastronomie Millimetrique (IRAM) Plateau de 



Bure Interferometer. The receivers were tuned to the para-H^O - 22,0 transition at 203.407498 GHz 
(1.47 mm). The correlator was set up with one unit with a bandwidth of 36 MHz (53 km s" 1 ) centered on 



3 In this paper we adopt a distance of 250 pc for NGC 1333-IRAS4B (see discussion in lEnoch et al. IJ200H)). 
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this frequency providing a spectral resolution on 460 channels of 0.078 MHz (0. 1 1 km s ). The source was 
observed in two configurations: in the C configuration on 02 December 2008 and in the B configuration 
on 11 and 13 January 2009. About 11 hours were spent in each configuration (including time used on gain 
calibrators etc.). When combined these two configurations cover baselines with lengths from 17.8 to 452 m 
(12 to 308 kf). 

The data were calibrated and imaged using the CLIC and MAPPING packages from the IRAM GILDAS 
software. The calibration followed the standard approach: the absolute flux calibration was established 
through observations of MWC 349, the bandpass by observations of the strong quasar 3c454.3 and the 
complex gains by regular observations of the nearby quasar J0336+323 (approximately 0.8 Jy at 1.45 mm). 
Integrations with clearly deviating amplitudes and/or phases were flagged and the continuum was subtracted 
prior to Fourier transformation of the line data. With natural weighting the resulting beam size is 0.67" x 
0.55"at a position angle of 36.7°; the field of view is 25" (HPBW) at 1.45 mm. The resulting RMS noise 
level is 11 mJy beam -1 channel -1 for the line data using natural weighting. The continuum sensitivity is 
limited by the dynamical range of the interferometer and the resulting RMS noise level is a few mJy beam -1 . 



3. Results 



Fig-fJJshows me continuum image of the observed region around IRAS4B. As seen both IRAS4B and its 
nearby companion IRAS4B' are clearly detected in the image. Table[T]lists the results of elli ptical Gaussian 
fits to t he two sources: both are resolved with fluxes in agreement with the results from IjQrgensen et all 
(|2007ah assuming that it has its origin in thermal dust continuum emission with F v oc y a with a « 2.5 - 3. 
The continuum peak is clearly offset by 5-7" from the emission at 3.6-24 fim seen in the Spitzer Space 
Telescope images of IRAS4B ; an indication that the Spitzer emissio n has its orig in in material hea ted by the 
protostellar outflow even at 24 /urn (see also ljprgensen et al.ll2007bl and Fig. 2 of [Alien et al.ll2007l) . 



Fig. [2] shows the spectrum toward the continuum peak of IRAS4B. A number of lines are clearly de- 
tected as listed in Table [3]- including the targeted H^O 3i % - 2?n line. For the line identification we used 
the JPL ^Pickett et al.ll 1998b and CDMS (|Muller et alJboOll 12005 ) databases and cross-checked those with 
the online Splatalogue compilation. Since all of the lines are narrow little line-blending occurs, in con- 
trast with high-mass star-forming regions. Also, all the assigned lines are centered within 0.1-0.2 km s -1 



Table 1: Parameters for IRAS4B and IRAS4B' from elliptical Gaussian fits to their continuum emission. 

IRAS4B IRAS4B' 
Flux 0.586 Jy 0.128 Jy 

R.A. (J2000) 03:29:12.01 03:29:12.84 
DEC(J2000) 31:13:08.07 31:13:06.93 
Extent" 0.80"x0.54" (-65°) 0.56"x0.45" (-86°) 

"Size of Gaussian from fit in (u, v)-plane (i.e., deconvolved FWHM size) and position angle of major axes 

(in parentheses). 



03 h 29 m 13!0 03 h 29 ,B 12!5 03 h 29 m 12?0 03 h 29 m ll!5 

a(J2000) 



Fig. 1. — Continuum image of IRAS4B and IRAS4B' from the IRAM PdBI observations (contours shown 
in logarithmic steps from 0.015 to 0.25 Jy be am 1 overlaid on Spitzer Space Telescope mid-infrared image 
(IJ0rgensen et alj|200q : iGutermuth et al.l I2008T) with 4.5 //m emission shown in blue, 8.0 /mi in green and 
24 /mi in red. 
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of the systemic velocity of 7.0 km s of IRAS4B. Most prominent in the spectrum are lines of dimethyl 
ether, CH3OCH3, with 5 identified transitions. In addition, lines of sulfur dioxide, SO2, and water, H,j, 8 0, 
are clear detections with a fainter line of ethyl cyanide, C2H5CN, also present. Concerning the confi- 
dence of the assignment of the H' 8 line: according to the Splatalogue compilation no other transitions fall 
within +1 km s" 1 of the location of the Hi, 8 line. Offset by 1-2 km s" 1 are transitions of 13 CH 2 CHCN, 
13 CH3CH2CN and (CH^^CO: the two former can be ruled out because of the lack of additional components 
which should have been observable at larger velocity offsets, whereas the latter has a very low intrinsic line 
strength. 

For each of the detected lines we fit a circular Gaussian in the (u, v) plane to the line emission integrated 
over the widths of the lines to zero intensity (about ±1 km s -1 ). The derived fluxes, peak positions and sizes 
are given in Table [3] as well. All the detected lines are marginally resolved with deconvolved sizes of 0.2"- 
0.6", corresponding to radii of only 25 to 75 AU. Interestingly, most of the detected lines are offset by about 
0.10-0.15"from the continuum position in right ascension - with the exception of C2H5CN. The peaks of 
the molecular emission still fall within the deconvolved extent of the continuum emission, though. The 
images of the complex organic molecules and SO2 will be discussed elsewhere. 



4. Discussion 



4.1. Velocity field 



The high spectral and spatial resolution offered by the IRAM data reveals a tentative velocity gradient 
in the H^O emission. Fig. [3] shows a moment-1 (velocity) map of this emission indicating a change in 
velocity of a few x 0.1 km s _1 over the extent of the source emission. To estimate the magnitude of the 
velocity gradient the centroid of the emission was determined channel-by-channel in the (u, v)-plane. The 
derived offsets along and perpendicular to the largest velocity gradient are then plotted as function of velocity 
and a linear fit performed (Fig. |3]). The velocity gradient found in this way is 9.4 km s _1 arcsec -1 or 
7.7xl0 3 km s" 1 pc~ ! at a 3.5^-cr confidence level. The velocity gradient is abo ut 3 orders of magnitu de 



larger than the typical gradients observed on arcminute scales in pre-stellar cores (IGoodman et al.lll993l) 



supporting a scenario in which the HJ, 8 emission has its origin in a more rapidly rotating structure such as 
a central disk - although it is not possible to address whether the velocity for example is Keplerian in nature. 
The position angle of the largest gradient is «75° measured East of North (with an accuracy of ±5-10° 
from the gradient measured in the (u, v) plane), which, interestingly, is nearly perpendicular to the axis of 
the protostellar outflow (upper panel of Fig. |3]) traced by water masers in the North-South direction with a 



position angle of -29° (Marvel et al 



at a position angle of 0° (e.g., 



Choi 



2008 



2001 



Desmurs et al. 



Jorgensen et al.ll2007ah . 



2009) and thermal line emission on larger scales 



Besides the tentative velocity gradient, the Hj, 8 line width is remarkably narrow: its width to zero 
intensity is about 1.7 km s" 1 with a full width half maximum (FWHM) from a Gaussian fit of 1.0 km s _1 . 
For comparison, a free-falling envelope toward a central star with even a very small mass of 0. 1 M Q would 
have a characteristic infall velocity of 2.5 km s -1 at 25 AU - i.e., we would expect a line width significantly 
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Table 2. Identified lines and results from Gaussian fit to emission in (u, v) plane. 



Molecule 


Transition 


Frequency 


E u 


/u 2 S 


Flux a 


Offset (Ace, A6) b 


Size c 






[GHz] 


[K] 




[Jy km s^ 1 ] 


["] 


["] 




3 1,3 - 22,0 


203.407498 


203.7 


0.344 


0.081 


(-0.11 


-0.02) 


0.20" 


CH3OCH3 


3 3 ,i - 2 2 ,i EA 


203.402779 


18.12 


9.29 


0.035 


(-0.14 


-0.07) 


0.48" 




33,i - 22.1 EE 


203.410112 


18.12 


31.1 


0.088 


(-0.14 


-0.02) 


0.46" 




3 3 ,o - 2 2 ,i AE 


203.411431 


18.12 


20.0 


0.044 


(-0.12 


-0.02) 


0.31" 




3 3 ,o - 2 2 ,i AA 


203.418656 


18.12 


33.4 


0.12 


(-0.12 


+0.03) 


0.67" 




33,o - 22,1 EE 


203.420253 


18.12 


22.3 


0.094 


(-0.16 


+0.00) 


0.25" 


C2H5CN 


232,22 - 222,21 


203.396581 


122.3 


338 


0.045 


(+0.02 


+0.00) 


0.37" 


S0 2 


12o,i2 - Hi, 11 


203.391550 d 


70.12 


22.5 


0.12 


(-0.11 


-0.02) 


0.41" 



a Total flux from fitting circular Gaussian to integrated line emission in the (u, v) plane. For conver- 
sion to brightness temperatures, the gain of the interferometric observations with the current beam size is 
0.0126 Jy K _1 - i.e., implying an integrated line strength for the H^O line of 6.5 K km s _1 . 

b Peak offset with respect to continuum peak estimated from Gaussian fit. 

c Extent of emission (FWHM) from Gaussian fit. 

d Catalog frequency uncertain (accuracy ±0.1MHz). The observed SO2 peak indicates a frequency higher 
by 0.1 MHz than the tabulated one. 
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Fig. 2. — Le/f.- Integrated emission of the Hi, O line (contours in steps of 5 mJy beam" km s starting at 
10 mJy beam" 1 km s -1 ) plotted over the continuum emission (grey-scale). Right: Spectrum extracted in the 
central 0.6" x 0.5" beam toward the continuum position for IRAS4B. The detected lines are indicated at 
the position of their catalog rest frequency corrected for the 7.0 km s" 1 systemic velocity of IRAS4B. The 
spectrum has been binned to twice the observed resolution The insert shows a blow-up of the Hi, 8 line at 
the original resolution. 
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6.0 6.5 7.0 7.5 8.0 

Velocity [km s -1 ] 




Fig. 3. — Velocity gradient in the Hi, O emission toward IRAS4B: Upper panel: Moment-1 (velocity) map 
with the direction of the largest gradient indicated by the solid line and the direction of the water masers 
with the dashed line and the larger-scale molecular outflow with the dotted line. Lower panels: the fitted 
position offset (relative to the phase center of the observations) along (left) and perpendicular (right) to the 
direction of the main gradient for each channel where at least 3cr emission is detected. The solid lines in 
both panels indicate the continuum position; the dashed line in the left-hand panel a linear gradient fit to the 
data. 
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larger than that observed in this case. A rotationally supported edge-on disk would show a similar line 
width, whereas a more face-on disk, would provide a lower width in agreement with what is observed, 
although a disk seen entirely face-on would not produce any obs ervable velocity gradient. Alternatively, 
an inc lined disk with sub-keplerian rotation (e.g., a "pseudo-disk" (iGalli & Shulll993l : lHennebelle & Ciardi 



2009)) or a densit y enhancement in the inner envelope due to a magnetic shock wall slowing down the 



infalling material (|Tassis & Mouschovias! 120051: IChiang et al.l 120081) ) could result in a l ower line width - 
although such phenomena usually occur on larger scales than 25 AU. IWatson et al.l (|2007l ) argued based on 
the detected mid-IR emission that the outflow cone must be close to pole-on providing a low opacity route 
for the IR emission to escape toward the observer. This would imply a nearly face-on disk. However, the 
prop er motions and rad i al velocities of H2O masers suggest an outflow directed closer to the plane of the 
sky ( Marvel et al. 2008 : Desmurs et al. 20091 ) . Clearly observations of the dynamical structure of the inner 
regions of IRAS4B with even higher angular resolution are required to distinguish between these scenarios. 



4.2. Column density, mass and abundance 



The strength of the Hi, 8 transition suggests a significant reservoir of water vapor toward IRAS4B. We 
estimate the colum n density of H2O by adopting an excitation temperature for p-Hi, 8 of 170 K from the 



modeling results of IWatson et all (120071) and assume that the emission is optically thin and uniform over its 



extent. With these assumptions we estimate a column density for p-H^O of 4 x 10 15 cm 2 , which translates 



into a total H2O column density of 9 x 10 cm assuming a 0/ O ratio of 560 and an ortho-para 
ratio for H2O of 3 (i.e., an ortho-para ratio establish ed at high tempera tures). For comparison the column 
density estimated based on the Spitzer detections by IWatson et al.1 (|2007l) is two orders of magnitude lower, 



9.2x 10 16 cm 2 , over an emitting area of 0.24"x0.24" - i.e., comparable to the deconvolved size of the H* 8 
emission here. Our inferred column density is almost unchanged if the excitation temperature is lowered to 
100 K and increases by up to a factor 5 if the temperature is increased to 1000 K. These temperatures cover 
the range of conditions expected for any of the scenarios for the origin of the H2O emission discussed in £l] 

The total mass contained in the detected H2O is 6.0xl0~ 8 M s , or 0.02 M^anh- Assuming a typical 
H?Q abundance relative to H2 of 10~ 4 , corresponding to sublimation of the H2O rich dust ice mantles (e.g., 
Pontoppidan et al.l 120041 ). the total H2 mass of the H2O emitting material (dust+gas) is 7.5xl0~ 5 M© or 
25 Mearth- For comparison the mass of the compact disk aroun d IRAS4B inferred fro m the modeling of 



high angular resolution dust continuum observations is 0.24 M© dJOrgensen et al.ll2009l ). Thus, if the HJ, 8 



emission has its origin in this disk, it arises in a small fraction « 0.03% of the material in the disk. 

Alternatively, in the absence of such a disk, it is possible that the emission has its origin in the hot inner 
region of the protostellar envelope where the temperature is > 100 K: for a simple power-law envelope den- 
sity profile reproducing th e submilhmeter contin uum emission for IRAS4B on scales larger than ~1000 AU 
(2.8 M© within 8000 AU: IjQrgensen et all 120091) . the mass within 25 AU (where the temperature is higher 
than about 100 K) is about 5xl0~ 4 M©, implying a H2O abundance of about 1.5xl0~ 5 . However, such 
a model is not self-consistent on small scales: to fit the observed compact dust continuum emission seen 
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by the interferometer a strong increase in the envelope density on smal l scales by two orders of magnitude 
is required - above the already increasing radial density profile (e.g., IjQrgensen et al.ll2009l). Howev er, if 
such a density enhancement was due to a magnetic field wall as discussed above (iChiang et al.l l2008). the 
H2O abundance would be lower by the same amount, dropping to about 1.5xl0~ 7 . This abundance is low 
compared to the expectation from the full desorption of the H 2 mantles and also lower than the constraints 
on the H2O a bundance in the o uter envelopes of the IRAS4 sources where H2O is frozen-out based on ISO- 
LWS results ( Maret et"al]l2002). A low H?0 abun dance in the region of grain-mantle desorption may reflect 
destruction of H2O by X-rays Jstauber et alJboOtih . but the H2O abundance would need to be reduced to the 
levels of the outer cold envelope where H 2 is frozen-out and thus could not provide the compact emission 
observed here. 

Models of the chemistry in more evolved disks around pre-main sequence stars (where the envelope 
has dissipated) show a warm upper layer where H2O gas can exist. Although these models are not fully 
appropriate for disks in the embedded phase, where UV photons may not be able to freely reach the disk 
surface and heat the gas, they provide a useful reference point for comparison. If the H2O gas phase abun- 
dance is just determined by the balance of photodesorption of H2O ice and freeze-out, typical gaseous H2O 
column densities at 10-25 AU are a few xlO 18 cm" 2 , dropping rapidly at larger radii in these models (e.g., 
Oberg et al.ll2009T) . only slightly lower than those found here. Alternatively, th e temperature s in the upper 
layers of the disk may be hot enough to drive an active gas-phase chemistry. IWoitke et al.l (12009b find a 
layer of irradiated hot water at altitudes z/R = 0.1-0.3 extending out to 30 AU where tempera tures are 200- 



1500 K and densities 10 8 - 10 10 cm 3 , comparable to the conditions deduced here (see also lGlassgold et al. 
(|2009|)). The H2O mass in this layer is ~ 10 A/earth in their model, about two orders of magnitude lower 



than derived on basis of the H^O observations presented in this paper. Because their H2O/H2 abundance is 
only 10~ 6 - 10~ 5 , the inferred total warm H2 mass is comparable. 

The disc repancy between the column densities from these data and those from the Spitzer mid-infrared 
observations (IWatson et all 120071) remains significant, though. A possible explanation is that the mid- 
infrared observations are limited by extinction and thus do not probe the total water column density. Alter- 
natively, a lower temperature of the H2O emitting gas of ~ 100 K could with an unchanged column density 
keep the observed Hi, 8 line intensity at the same level while decreasing the mid-IR line flux predicted by 
the model. 



It a lso remains a quest i on why IRAS 4B is the only source with strong H2O lines at mid-infrared wave- 



lengths (IWatson et al. 2007). 



Watson et al. 



proposed that IRAS4B was either seen in a particularly favorable 
geometry - with the H2O mid-infrared emission from disks around other sources being masked by optically 
thick envelopes - or that the release of H2O into the gas-phase constituted a particularly short-lived stage 
in the evolution of embedded protostars. These scenarios could easily be distinguished by further observa- 
tions of H^O emission from embedded protostars at (sub)millimeter wavelengths where the envelopes are 
optically thin and water should thus be detectable if present at similar levels as in IRAS4B. Future obser- 
vations with the Atacama Large Millimeter Array in Band 5 will fully open this topic up for ground-based 
observations of large samples of embedded protostars in the same lines. 
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